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ABSTRACT 

For 18 well-observed gravitationally lensed QSOs, we compare new non-parametric mass profiles for the lens- 
ing galaxies with stellar-population models derived from published HST photometry. The large volume of param¬ 
eter space searched - with respect to the possible star formation histories - allows us to infer robust estimates and 
uncertainties for the stellar masses. The most interesting results are: (1) the transition from little or no dark matter 
in the inner regions (<j r e ) to dark matter dominating on the ~ 5 r e scale (~ 20kpc) is clearly seen in massive el¬ 
lipticals; (2) Such a trend is not seen in lower-mass galaxies, so that the stellar content dominates the mass budget 
out to ~ 5 r e \ (3) the radial gradient in the dark-matter fraction for these intermediate redshift galaxies agrees with 
published data on nearby galaxies. This result can help reconcile the discrepancies found in recent estimates of 
dark matter in elliptical galaxies using different techniques (e.g. Planetary nebulae versus X-ray). The observed 
trend suggests the stellar component in massive galaxies extends further out in terms of the dark matter scale 
radius. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: elliptical and lenticular, cD — dark 
matter 


1. INTRODUCTION 

For galaxies further than ~ lOOMpc, the usual kinematic 
tracers of galactic dynamics become increasingly difficult or 
impossible to use. But at larger distances, nature sometimes 
provides a very different indicator of galaxy mass—strong lens- 
ing of quasars. Galaxies with a quasar conveniently placed be¬ 
hind them, which they then lens into multiple images, are rare. 
But for those galaxies it is fairly easy to measure masses, even 
to z ~ 1. In recent years there has been progress on estimat¬ 
ing the M/L in samples of lensing galaxies (Keeton et al. 1998, 
Kochanek et al. 2000, Rusin et al. 2003). For a few systems 
there have been efforts to combine lensing and velocity disper¬ 
sions (Treu & Koopmans 2002a,b, Koopmans & Treu 2003). 

In this Letter we go beyond simple M/L for a sample of lens¬ 
ing galaxies and try to recover the distribution of stellar and 
dark mass within galaxies. Our technical innovations are (i) we 
use observed colors to model the stellar population in detail and 
hence map the stellar mass, and (ii) we use the method of pix- 
elated lens reconstruction to make detailed profiles of the total 
mass; we pay particular attention to quantifying the uncertain¬ 
ties in both departments. 

Our sample comprises 17 early-type galaxies over a wide 
range of redshifts (0.3 < z < 1) and a bulge (z = 0.04). The 
sample has been selected from the CASTLeS group database 1 . 
The main properties are listed in Table 1. We assume a concor¬ 
dance cosmology (O m = 0.3, Oa = 0.7) with H / 1 = 14Gyr. 

2. MEASURING LENSING MASSES 

For most lensed quasars, it is fairly easy to fit a model galaxy 
lens to the image positions and hence provide a map of the total 
sky-projected mass. But such a map, although a reasonable first 
approximation, will be non-unique because of lensing degen¬ 
eracies (Falco et al. 1985, Gorenstein et al. 1988, Saha 2000). 
A way around this problem is to generate an ensemble of mod¬ 


els (Williams & Saha 2000, Trotter et al. 2000, Keeton & Winn 
2003) all models constrained to reproduce the observed image 
positions (and time delays if known) precisely. From the model 
ensemble, estimates and uncertainties of any desired quantity— 
for example, the mass at a given projected radius—can be easily 
extracted. We use the PixeLens code (Saha & Williams 2004) 
which more or less automates the whole procedure, even for 
complex lenses like B1608+656. 

The model-ensemble technique (we used 200 models per 
lens) makes our mass uncertainties larger than in previous work 
on M/L ratios, but much more realistic. Four-image lenses 
tend to be better constrained than two-image lenses. Also, sys¬ 
tems with known time delays allow much tighter mass estimates 
(provided Hq is assumed, as we do here). 

3. MEASURING STELLAR MASSES 

The stellar mass content can be determined from the photom¬ 
etry although a fair share of assumptions must be invoked in or¬ 
der to transform light into mass. These assumptions relate to the 
age, metallicity and mass distribution of the unresolved stellar 
populations. The latter can be reduced to a time-independent, 
universal initial mass function (IMF) as suggested by obser¬ 
vations (see e.g. Wyse 1998). In this paper we explore two 
IMFs: Salpeter (1955) and Chabrier (2003), both defined be¬ 
tween O.IMq and IOOMq. The former assumes a simple power 
law over the allowed range of stellar masses. The IMF proposed 
by Chabrier (2003) has a very similar upper-mass dependence 
as the Salpeter function. However, the low-mass end assumes a 
flatter - more physical - behaviour, following a lognormal dis¬ 
tribution. This IMF gives M/L ratios which are a factor ~ 1.5 
smaller than those for a Salpeter IMF. 

The main difference between most of the IMFs proposed 
(e.g. Salpeter 1955; Miller & Scalo 1979; Kroupa, Tout & 
Gilmore 1993; Chabrier 2003) lies in the relative contribution 
from stars with masses M < M 0 . Since these stars have a large 


See the CASTLeS page, http: //cfa-www. harvard. edu/glensdata/ for a list. 

1 


2 




012345 01234 

r / r e r / r e 

Fig. 1 . — Profiles of the mass enclosed as a function of radius for two low-mass (left) and two massive early-type galaxies (right). 
The solid and empty circles give the total and stellar mass content, respectively. The stellar masses assume a Chabrier IMF but a 
Salpeter IMF leads to qualitatively the same result. The error bars indicate 90% confidence intervals. The vertical dashed line gives 
the position of ri ens . 


M/L ratio, stellar populations born with either IMF feature very 
similar spectral energy distributions, albeit with a different nor¬ 
malization for a given stellar mass (see e.g., Bruzual & Char¬ 
iot 2003). Therefore, a photo-spectroscopic analysis cannot be 
used reliably to constrain the IMF at the low-mass end. In¬ 
stead, we will use the available photometry from our sample 
and compute the stellar mass content corresponding to either a 
Salpeter or a Chabrier IMF. These two mass distributions repre¬ 
sent a robust range of possible values between the unphysically 
high M/L imposed by a simple power law extrapolated to low 
masses and a more realistic distribution. 

The age and metallicity distribution can be constrained by 
a comparison of the photometry with a simple model of star 
formation. We assume a 3-parameter model which reduces 
the description of the stellar populations to a single metallic¬ 
ity (Z*), a formation epoch (/forX and a formation timescale 
(tsf), so that at any given time the star formation rate is rp(t) oc 
exp(-At /tsf), with At = f-fFOR- Each choice of parame¬ 
ters (Z*, ?for i tsf) represents a possible formation scenario, and 
can be convolved with simple stellar populations (SSP) in or¬ 
der to generate a composite model from which various photo- 
spectroscopic observables can be retrieved. We use the Bruzual 
& Chariot (2003) population synthesis models. 

We use the published photometry of the lenses from the 
CASTLeS group (Rusin et al. 2003) which was obtained from 
HST/WFPC2+NICMOS images after carefully subtracting the 
contribution from the sources. Each set of measurements cor¬ 
respond to a number of colours (from 1 to 5) which is used in 
order to constrain the parameter space described above. We use 
their radial fits to generate the profile of the stellar component. 
Galactic reddening was included in the analysis by using the 
E(B-V) values from Rusin et al. (2003) and applying a dust 
correction according to the R=3.1 curve of Fitzpatrick (1999). 

The best fit for each lens was obtained using an adaptive grid 
in the 3-dimensional parameter space that describes all possible 
star formation histories in our model, and using a Metropolis al¬ 
gorithm (see e.g., Binney et al. 1992, or Saha 2003) to find the 


uncertainties which are quoted with respect to the 5th-95th per¬ 
centiles (90% confidence interval) throughout the paper. The 
stellar mass content is computed for the best fit, with respect to 
the light inside the half-light radius ( r e ). 


4. STELLAR VS TOTAL MASS 

Figure 1 shows the mass profiles (stellar and total) for four 
galaxies, two with M tot 10 I2 Mq (left) and two with M tot :> 
10 i2 Mq (right). The result is striking: the low-mass galaxies 
have little or no dark matter at all observed radii. The high- 
mass galaxies have little or no dark matter inside of r ~ r e , but 
at large radii they are dominated by dark matter. Here we have 
a clear indication that massive galaxies have a gradient in their 
dark-matter fraction. 

Rather than presenting mass profiles for all the galaxies in 
our sample, we will summarize the profiles in two ways: first 
by comparing stellar and total mass at a chosen radius, and then 
by quantifying the dark-matter gradient. 

In Figure 2 and Table 1 we show three derived aperture 
masses: the total mass derived from lensing, the stellar mass 
derived assuming a Salpeter IMF, and the stellar mass derived 
assuming a Chabrier IMF. Uncertainties are at 90% confidence. 
The aperture radius in all cases is a quantity we call ri em , 
and it is the largest radius at which the lensing mass is well- 
constrained. In general, ri ens represents the radius over which 
there is lensing information, and is roughly the radius of the 
outermost image. From Figure 2 it is apparent that the dark- 
matter fraction increases with mass, and the trend is compatible 
with the observed tilt of the NIR fundamental plane (Mobasher 
et al. 1999; Ferreras & Silk 2000). Additionally, we see that M st 
is estimated as more than M t0 , in a few cases for the Salpeter 
IMF, but never for the Chabrier IMF; while it would be cavalier 
to argue against a Salpeter IMF on this basis, it does suggest 
that future comparisons of lensing and stellar-population mod¬ 
els may help constrain the IMF. 

We now consider the gradient in the dark-matter fraction. 
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Fig. 2.— A comparison of total and stellar masses for our sam¬ 
ple of gravitational lens early-type galaxies (all measured inside 
rims). The upper and lower panels correspond to a Salpeter and 
a Chabrier IMF, respectively. The solid line represents M tot /M st 
and the dashed line follows the expected correlation from the tilt 
of the fundamental plane for ellipticals. 


Following Napolitano et al. (2005; hereafter N05) we introduce 



Even though VT is a noisy quantity, it has advantages over 
M, ot /M sl . Being a gradient it is not sensitive to the choice of 
wavebands used to measure stellar luminosity. Also, VT is 
unlikely to be very sensitive to the redshift evolution of galax¬ 
ies, because to first order M to , /M st will change at the same rate 
throughout the inner parts of a galaxy as its stellar population 
ages. We are able to compute VT with tolerable uncertainty 
for only a subsample of galaxies. 

The work of N05 differs from ours in several respects: 
(i) they derived masses from the kinematics of planetary nebu¬ 
lae and globular clusters; (ii) their galaxies are nearby whereas 
ours span range from z — 0.3 to and z — 1; (iii) they used three- 
dimensional radial gradients while we are working with pro¬ 
jected quantities. Point (ii) is unlikely to affect VT much, as 
we have already argued. Regarding point (iii), the difference is 
minimal for the estimation of the stellar masses, since as Fig¬ 
ure 1 shows, the stellar mass is concentrated within the central 
few kpc; even for the more extended distribution of the dark 
matter, most of the mass that makes up the total line of sight 
mass resides at small radii. In general, 2D and 3D VT esti¬ 
mates differ by less than 50%. We conclude that we can directly 
compare with N05. 

Figure 3 compares our estimates of VT with those of N05 
in the format of their Fig. 7. There is a general agreement 
between the two sets of results, at least over the mass range 
covered by both data sets. Our values are somewhat higher 
for 11 < log(M*/M 0 ) <11.5 range, but with the current small 


Fig. 3.— Dark matter gradient VT, computed between r,„ ~ 
0.25 r e and r out ~ 4.5r,,. Filled circles correspond to the Chabrier 
IMF and empty circles to the Salpeter IMF. The error bars cor¬ 
respond roughly to 1 cr. This figure may be compared with Fig. 7 
of Napolitano et al. (2005). The star symbols are taken from that 
paper. The curves are la bands for model predictions with dif¬ 
ferent star-formation efficiencies, also adapted from Napolitano 
et al. 

numbers we cannot assess significance. In Figure 3 we have 
also plotted the predictions of ACDM models by N05. This re¬ 
quires some justification. First, the models were computed for 
the r e -M st relation characteristic of the N05 data set, but our r e - 
M st turns out to be very similar. Second, the models are based 
on NFW ACDM predictions for halos at z ~ 0; do they still 
apply to our data, with z ~ 0.3- 1? Fortunately, yes, because 
the relation between the concentration parameter and halo mass 
is fairly robust with respect to the redshift of observation in 
the z = 0-3 range (see Fig. 14 of Wechsler et al. 2002). The 
curves in Figure 3 show three different star-formation efficien¬ 
cies: 6sf = 0.07,0.25, and 1, where e$F is defined as the ratio be¬ 
tween total stellar and baryonic matter. The general trend is that 
higher stellar mass systems should have a monotonically larger 
contribution from dark matter with increasing radius. This can 
be understood intuitively as follows (cf. section 3.4 of N05). In 
galaxies with higher stellar masses the stellar component ex¬ 
tends to larger radii, in terms of NFW profile scale radius, and 
so larger radii are more dark matter dominated. Our data is con¬ 
sistent with that trend. However, one should be cautious when 
interpreting the results with respect to 6sf ■ The N05 models 
disregard the effect of the gas collapsing to the centre of the 
halo and its effect on the dark matter density profile. 

5. DISCUSSION 

In the Milky Way it is generally held that in the inner few kpc 
and on mass scales of ~ 10 1 'M© there is little or no dark matter, 
whereas on the scale of ~ lOOkpc and ~ 10 12 M© dark matter 
completely dominates. In this Letter, we see the same transi¬ 
tion in elliptical galaxies at redshifts of 0.3 to 1, based on pro¬ 
jected mass distribution derived from lensing. The dark matter 
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Table 1 

Aperture Mass (x 10 10 M q ) 


Object 

z 

r f /kpc 

r /f „,,/kpc 

M v 

flens) 

M va /p(<C Vlens) 

M chab (^^7 ens ) 

B0047-280 

0.49 

5.4 

8.6 

-22.5 

70 '736.2 
'32.3 

n.ift 5 

9-Oft 7 

B0142-100 

0.49 

3.1 

5.6 

-22.7 

24.9 3 ‘; 7 

18.3ft! 

20.9??;? 

J0414+053 

0.96 

14.5 

29.0 

-23.1 

82.5 ? 3 ;? 

74 7115.4 
,Qc - ' 44.5 

^ 1 oso.s 
J1 * J 31.0 

B0818+122 

0.39 

4.8 

10.6 

-22.0 

67.4”;? 

20.8?3 4 

16.2?!;? 

J0951+263 

0.20 

0.7 

2.1 

-19.3 

4 7 5-7 
'3.6 

1-5?;? 

Lift 

B0952-011 

0.38 

0.5 

2.5 

-19.4 

4.5ft 

4.4ft 

0 r4.0 
J.3 2 .7 

B1009-025 

0.78 

1.6 

4.8 

-21.0 

18.0“;? 

749.8 

7 -^5.0 

5.5ft 

J1017-204 

0.86 

2.4 

2.4 

-21.6 

4.8“ 

6.4ft 0 

4.3ft 0 

B1030+074 

0.60 

2.6 

4.7 

-21.8 

16 . 1 ft? 

14.5ft 3 

10.6ft 3 

B1104-181 

0.73 

4.4 

16.3 

-22.7 

122.0”?;° 

36.6“;? 

22.8?!;? 

Bll 15+080 

0.31 

2.3 

7.4 

-21.1 

?7 S 31 - 5 
z/ - :) 25.4 

14.8ift 

10.5ft 3 

J1411+521 

0.46 

2.9 

9.3 

-21.2 

oo zr 38.8 
•^•^•^ 29.6 

10.0‘ft 

8.9ft 8 

J1417+522 

0.81 

5.2 

16.1 

-23.4 

118-Ojfto 

43.4ft? 

32.8??'? 

B1422+231 

0.34 

1.5 

6.0 

-20.7 

15-8 ft? 

7.9ft 3 

r ^)10.6 
6-2J.! 

B1520+530 

0.72 

2.2 

6.6 

-22.4 

26.9”;? 

21.8 3 ?;' 

18.5?°;? 

B1608+656 

0.63 

4.2 

6.3 

-23.1 

85.0?°;? 

40.2ft 3 

28.5? 7 ; 8 

B2149-274 

0.50 

2.9 

4.1 

-21.9 

13.9!?;? 

6.9?;? 

4.6ft 

B2237+030 

0.04 

3.1 

0.9 

-20.9 

L9?j 

L6ft 

L9ft 


fraction tends to increase with the total mass. Incidentally we 
note that one of the possible mechanisms that can be invoked to 
explain the tilt of the fundamental plane involves such a trend 
between stellar and total mass with a very similar scaling be¬ 
havior (Ferreras & Silk 2000). Part of the tilt is explained by 
non-homology effects (Trujillo et al. 2004), but our data sug¬ 
gest a further contribution to the tilt from the different stellar 
vs dark matter distributions. We claim that most of the tilt of 
the fundamental plane - at least when the observables used to 
define it reach beyond r e - is caused by the correlation between 
dark matter fraction and galaxy mass as presented in this Letter. 
The gradient of the dark matter fraction inside individual galax¬ 
ies is also similar to what is inferred in nearby galaxies. These 
results indicate that the distribution of dark matter in galaxies 
has not changed much since z = 1 • 

Enlarging the sample of galaxies for which stellar and dark- 


matter profiles can be reconstructed in this way is an obvious 
goal for future work. But another line of future work is to im¬ 
prove the spatial resolution in some cases. A possible outcome, 
hinted at in this work, is that bottom-heavy IMFs such as the 
Salpeter function might be ruled out by the comparison between 
stellar and lensing masses. 

The correlation presented in this Letter help to reconcile the 
apparent discrepancy between the low dark matter content from 
dynamical estimates using planetary nebulae (Romanowsky et 
al. 2003) 2 and the dark matter dominated galaxies suggested by 
X-ray data (Loewenstein & Mushotzky 2002). Mass estimates 
of faint X-ray ellipticals, albeit challenging, will definitively 
test the observed correlation. 

We would like to thank Nicola Napolitano for a very useful 
discussion and for providing his data and LCDM predictions. 
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